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Capillary  viscometer  measurements  were  maor  on  a series  of  concentrated 
suspensions  with  the  sane  fiber  length  d i st . i but  ion . Shear  stress-shear  rate 
data  are  reported  for  0,  15.  -<nd  3?  vol.?  glass  fiber  suspensions  in  100,  600,  and 
1000  P silicone  oil.  Yield  stresses,  wall  effects,  concentration  defects,  squeeze- 
through,  and  capillary  entrance  exclusion  were  not  significant  in  the  measurements. 

With  other  parameters  held  constant,  the  suspension  viscosity  increased  with  oil 
viscosity.  The  data  were  fit  by  a dimensionless  plot  of  the  form  F(fRe,De)  * 0 
where  f is  the  friction  factor  for  tube  flow  and  Re  and  De  are  the  Reynolds  and 
Deborah  numbers  of  the  system.  The  suspension  viscosity  increased  with  fiber 
concentration,  but  the  magnitude  of  the  increase  was  highly  shear  rate  dependent. 
The  data  are  compared  with  literature  predictions  for  the  dependence  of  suspension 
viscosity  on  fiber  concentration  and  shear  rate. 
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ABSTRACT 


Capillary  viscometer  measurements  were  made  on  a series  cf  concentrated 
suspensions  with  the  same  fiber  length  distribution.  Shear  stress-shear 
rate  data  are  reported  for  0,  15,  and  30  vol.i  glass  fiber  suspensions  In 
100,  600,  and  1000  P silicone  oil.  Yield  stresses,  wall  effects,  con- 
centration defects,  squeeze-through,  and  capillary  entrance  exclusion  were 
not  significant  in  the  measurements. 

With  other  parameters  held  constant,  the  suspension  viscosity 
increased  with  oil  viscosity.  The  data  were  fit  by  a dimensionless  plot 
of  the  form 


F(f Re , De)  = 0 

where  f is  the  friction  factor  for  tube  flow  and  Re  and  De  are  the 
Reynolds  and  Deborah  numbers  of  the  system.  The  suspension  viscosity 
increased  with  fiber  concentration,  but  the  magnitude  of  the  increase 
was  highly  shear  rate  dependent.  The  da»a  are  compared  with  literature 
predictions  for  the  dependence  of  suspension  viscosity  on  fiber  concentration 
and  shear  rate. 
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RHEOLOGY  OF  CONCENTRATED  SUSPENSIONS  OF  FIBERS  IN  TUBE  FLOW: 

III.  SUSPENSIONS  WITH  THE  SAME  FIBER  LENGTH  DISTRIBUTION 

I.  INTRODUCTION 

Part  I (1)  of  this  series  reviewed  the  sparse  experimental  literature 
on  the  rheological  properties  of  concentrated  fiber  suspensions  in  tube 
flow.  Little  agreement  obout  the  rheological  properties  was  found  among 
various  researchers.  Part  II  (2)  reported  the  results  of  an  exploratory 
study  of  concentrated  fiber  suspensions  In  which  a strong  dependence  of 
rheology  on  fiber  length  distribution  was  observed.  The  separate  effects 
of  fiber  concentration,  oil  viscosity,  and  fiber  length  distribution  on 
suspension  viscosity  were  difficult  to  measure  because  these  variables  are 
compounded  during  mixing.  This  paper  describes  the  results  of  viscosity 
measurements  on  a series  of  suspensions,  all  of  which  had  the  same  fiber 
length  distribution. 


1 1 ■ VISCOSITY  THEORIES 

In  order  to  make  comparisons  with  our  data,  we  review  briefly  the 

theoretical  studies  which  provide  explicit  expressions  for  viscosity 

applicably  to  concentrated  fiber  suspensions  in  tube  flow.  Brodnyan  (3) 

predicted : 


n * exp 

r 


( 2.54  + 0.399  (rp-l ) 1 * * * * • /*8  4 
| l-k4 


} 


(1) 


In  this  equation  nr  is  the  reduced  viscosity,  4 is  the  volume  fraction, 
rp  is  the  length-to-diameter  ratio  of  the  fibers,  and  k is  a crowding  factor. 
The  numbers  0,399  and  I.A8  are  constants  which  were  evaluated  empirically. 
Even  for  the  data  Brodnyan  used  to  evaluate  the  numerical  constants,  the 
equation  will  fit  the  data  only  up  to  0.5  vol.$. 
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More  recently,  Ziegel  (4)  derived: 


Hill  ♦ 1 - Z 


S'.  W 


♦ 1 


(2) 


where  u is  on  interaction  or  cioseness-of -approach  parameter,  B is  a rate 
constant  for  th»  equilibrium  between  free  particles  and  agglomerates,  and 
Z is  he  degree  of  agglomeration.  Only  the  shear  rate  dependence  of  the 
equation  was  tested  in  Ziegel *s  paper  and  the  test  was  only  over  a shear 

rate  range  of  one  and  a half  oroers  of  magnitude. 

[.  Hushir  (SI  used  the  flow-elasticity  analogy  to  predict: 
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This  evasion  is  expected  to  be  valid  only  for  parallel,  randomly  placed, 
infinitely  Ion,  fibers.  Nielsen  (6)  used  the  same  analogy,  with  the 
Hal p i n-Tsa i equation,  to  predict: 

1 ♦ (k-i)B4  (4) 
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where  <$m  is  the  maximum  possible  volume  fraction.  The  quantity  B Is  very 
nearly  unity  :or  rigid  particles  in  a viscous  medium.  To  the  authors' 
knowledge,  the  equation  has  never  been  tested  for  concentrated  fiber 
suspensions. 

Ill  . EXPERIMENTAL  TECHNIQUES 

The  extruder  mixing  technique  and  the  capillary  viscometer  used  for 
viscosity  measurements  in  this  study  were  described  in  detail  elsewhere  (2). 
The  samo  fiber  length  distribution  was  produced  in  all  suspensions  by 
uti  izing  the  fact  that,  with  additional  mixing,  the  suspensions  approached 
an  "ultimate  state  of  mixedness"  beyond  which  continued  mixing  at  the 
same  conditions  did  not  change  the  suspensions'  properties.  The  suspensions 
used  in  this  study  were  mixed  as  follows.  The  viscosities  of  all  silicone 
ous  were  nearly  equalized  by  changing  the  temperature  of  the  extruder. 

This  step  ensured  that  all  suspensions  with  the  same  fiber  concentration 
would  have  the  same  fiber  length  distribution  after  mixing.  The  suspensions 
of  highest  concentration  were  then  mixed  to  their  "ultimate  state  of 
mixedness"  (six  passes  through  the  extruder  were  required).  Suspensions 
of  lower  concentration  were  made  without  changing  the  fiber  length 
distribution  by  blending  the  high  concentration  suspensions  with  additional 
oil  and  extruder  mixing. 

The  mixing  procedure  afforded  little  independent  control  over  the 
fiber  length  distribution  produced;  but,  as  the  log-normal  distribution 
statistics  in  Table  1 show , the  distribution  was  nearly  the  same  in  all 
suspensions.  Figure  I shows  a typical  fiber  histogram;  the  mean  fiber 
length  is  0.005"  and  the  median  length  is  0.01". 


I 

! 

I 

i i 

t 

L* 

r 

r ... 

l 

» 


; 


L 


i 


1 


! 


I 


I 


b 

I V . SUSPENSION  PARAMETER 

We  believe  Table  2 lists  all  the  suspension  parameters  (as  opposed 
to  flow  parameters)  which  determine  the  viscosity  of  a concentrated 
suspension  of  fibers.  The  values  of  fiber  diameter,  modulus,  and  Poisson 
ratio  listed  in  Table  2 are  manufacturers1  specifications.  Since  these 
values  are  common  to  all  our  suspensions,  their  independent  determination 
was  felt  to  be  unnecessary. 

The  values  of  the  remaining  parameters  were  measured  as  follows.  The 
suspensions  were  mixed  by  weight,  so  the  fiber  concentration  was  determined 
at  the  time  of  mixing.  The  actual  concentration  achieved  was  verified 
for  selected  suspensions  by  Schxlet  extraction  of  the  oil.  Compressibility 
measurements  indicated  the  amount  of  residual  trapped  air  (the  third  phase 
of  major  concern)  was  less  than  0.12  vol.fc--a  negligible  amount.  In 
numerous  microscopic  examinations,  the  silicone  oil  was  observed  always 
to  wet  the  glass  fibers  and  the  fibers  were  observed  always  to  be 
separated  from  their  original  bundled  configuration.  These  observations 
indicate  negligible  surface  effects  and  random  disperson  of  fibers.  The 
fiber  orientation  is  believed  to  be  random  as  a result  of  the  high  shear 
in  the  extruder. 

The  oil  viscosities  as  functions  of  shear  rate  are  shown  In 
Figures  5~7 • The  oils  show  non-Newtonian  behavior  at  high  shear  rates. 
Possible  cahnges  in  the  oil  viscosity  during  mixing,  due  either  to 
degradation  or  to  evaporation  of  the  low  molecular  weight  fraction,  were 
monitored  by  centrifuging  the  oils  from  the  suspensions  and  measuring 
their  viscosities.  Although  the  viscosities  did  appear  to  change  a few 
percent,  experimental  difficulties  in  removing  ail  the  qlass  particles 
made  the  results  ambiguous  (7).  The  suspected  changes  were  small  and  no 
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correction  was  made  for  them  in  our  data. 
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V.  EXPERIMENTAL  RESULTS 

As  shown  in  Table  3,  no  differences  in  suspension  viscosity  or 
entrance  pressure  losses  were  observed  when  the  entrance  angle  of  the 
viscometer  was  changed  from  l80°  to  90°.  This  result  has  been  observed 
before  (8-10)  and  is  presumably  due  to  the  formation  of  a stagnant 
region  in  the  180°  entrance  such  that  the  effective  angle  is  closer  to 
90°.  All  other  measurements  reported  in  this  paper  were  made  with  the 
90°  entrance  angle. 

Shear  stress-apparent  shear  rate  data  for  the  suspensions  are 
plotted  in  Figures  2-4.  The  symbols  in  these  figures  represent  + 6* 
error  flags,  which  constitute  one  standard  deviation  in  our  measurements. 

None  of  the  suspensions  exhibited  a yield  stress  in  the  shear  rate  range 
measured  as  indicated  by  the  absence  of  a horizontal  asymptote  at  low 
shear  rates.  As  muy  be  seen  in  Figures  2-4,  no  differences  between  the 
flow  curves  measured  in  the  1/4"  and  1/8"  diameter  capillaries  were 
observed.  This  observation  precludes  the  possibility  of  measurable  wall 
effects,  such  as  si ip-at-the-wall , mechanical  exclusion,  or  hydrodynamic 
interaction;  probably  due  to  the  small  fiber  lengths  relative  to  the  tube 
diameters.  The  observation  makes  apparent  viscosity  a useful  parameter  for 
characterizing  the  flow  of  suspensions  in  tubes  with  diameters  greater 
than  1/8".  Apparent  viscosity-apparent  shear  rate  data  are  plotted 
in  Fiaures  5”7. 

Capillary  entrance  exclusion  and  squeeze-through  were  not  observed 
to  occur  wi  h any  or  the  suspensions.  The  possible  existence  of  a concen- 
tration defect  was  checked  by  comparing  the  capillary  and  ex.t  concentrations 
of  a suspension  in  the  flow  condition  mo'jt  likely  to  reveal  a defect: 


highest  flow  rate  and  longest  capillary.  No  concentration  defect  was 
observed.  These  results  are  not  surprising  in  view  of  the  absence  of 
measurable  wall  effects. 

VI . DIMENSIONAL  ANALYSIS 

In  order  to  discuss  the  effect  of  oil  viscosity  and  fiber  concentration 
on  suspension  viscosity,  it  is  useful  to  apply  dimensional  analysis  to 
our  system.  Dimensional  analysis  is  possible  with  the  aid  of  one 
assumption  whose  plausibility  we  shall  now  motivate.  We  restrict  our 
analysis  to  suspensions  with  the  same  fiber  length  distribution  and 
fiber  concentration  and  to  measurements  made  with  the  same  viscometer. 

We  treat  the  suspension  as  a continuous,  two  phase  system  of  rigid  fibers 
in  a generalized  Newtonian  fluid  and  assume  that  only  hydrodynamic 
forces  operate.  The  system  we  wish  to  non-d imensiona 1 i ze  includes  all 
of  the  suspension  in  the  viscometer  from  the  time  the  piston  starts  to 
move . 

The  major  difficulty  in  dealing  with  this  system  is  our  lack  of 
knowledge  about  the  initial  positions  of  the  fibers  where  the  boundary 
conditions  for  the  problem  are  given.  We  have  observed,  however,  that  if 
different  batches  of  the  same  suspension  are  run  through  the  viscometer, 
the  same  force  is  always  required  to  extrude  at  a given  piston  speed.  (2) 
Since  exact  duplication  of  the  initial  fiber  positions  from  batch  to 
batch  is  highly  Improbable,  we  conclude  that  there  must  be  many  different 
initial  distributions  which  are  statistically  equivalent  as  far  as  the 
average,  macroscopic  properties  measured  are  concerned.  This  conclusion 
seems  intuitively  reasonable,  for  the  initial  fiber  positions  arc  presumed 
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random,  and  there  are  many  possible  random  d i s tr i but  Ions . Furthermore, 
the  distributions  which  are  possible  do  not  depend  on  the  oil  viscosity, 
so  it  seems  reasonable  that  the  distributions  which  are  actually  produced 
during  mixing  also  may  not  depend  on  the  oil  viscosity  (over  a reasonable 
range  of  viscosities).  We  therefore  assume  that  two  well  mixed  suspensions 
with  the  same  f i be  length  distributions  and  fiber  concentrations  Put 
different  oil  viscosities  have  eqiival^nt  initial  fiber  position  distributions. 
Equivalent  in  this  sense  means  that  the  measured  viscosities  of  the  two 
suspensions  would  he  experimentally  i nd i s t i ngu • shabl e if  the  suspending 
fluids  were  the  same. 

With  the  above  assumption,  all  of  the  suspensions  under  consideration 
are  geometrically  similar.  Himensional  analysis  is  possible  since  the 
dimensionless  groups  which  charccerize  all  the  systems  are  the  same. 

To  determine  the  relevant  dimensionless  groups,  we  explicitly  non- 
d i mens iona I i ze  the  equations  of  motion  which  dascribe  our  suspensions. 

Although  any  t hree -constant  generalized  Newtonian  fluid  model  would  be 
acceptable  for  the  analysis,  we  choose  for  concreteness  the  Carreau 
mode  I Ml): 


i l r 


i ♦ (y/y0) 


n~rT-l 

J 2 


(5) 


The  Carreau  model  constants  for  the  silicone  oils  used  in  this  study  are 
listed  in  Table  4.  We  assume  that  the  constitutive  equation  which 
reptesents  the  fluid  in  viscorretric  flow  will  also  serve  to  represent  it 
in  whatever  flow  occurs  in  the  suspensions.  Elastic  effects  in  the 
silicone  oils  are  negligibly  small. 
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To  non-dimensional ize  the  equations  of  motion,  we  introduce  the 
following  dimensionless  variables: 

y*  - y/V  P*  « p4p-  l*  "T  “ x/D  (6) 

L o 

where  D is  the  capillary  diameter,  is  t»  e average  pressure  drop 

across  the  capillary,  and  V is  the  average  flow  rate  through  the  capillary. 

1 

The  constitutive  equation  may  be  written: 


o “ 


**)  i* 


(7) 


where  h is  a dimensionless  function  specified  by  the  Carreau  model. 
Neglecting  both  gravity  and  inertia,  the  equation  of  motion  becomes: 


0 = -EuRe  7*P*  ♦ 7*  • h(De,  n,  V)  Y* 

p -P  y 

The  dimensionless  groups  Eu  = — ~ , Re  * £ and  De  * jjy  are 

pV  o o 

the  Euler,  Reynolds,  and  Deborah  numbers,  respectively,  with  their 
attendant  physical  meaning. 

As  shown  in  Table  U the  power  law  index,  "n",  is  nearly  constant 
for  the  j i i i cone  oils,  so  we  ignore  explicit  dependence  on  "n". 

Dimensionless  qroups  representing  fiber  concentration,  diameter,  and 
length  distribution  are  not  explicitly  required  because  of  the  assumption 
that  these  quantities  are  constant  for  the  suspensions  considered.  One 
additional  dimensionless  group,  L/D,  the  ratio  of  length  to  diameter  of 
the  capillary,  is  necessary  to  meet  the  boundary  conditions.  The  suspensions 
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are  thus  characterized  by  the  three  dimensionless  groups  EuRe,  De,  and 
L/D.  Since  the  dimensionless  groups  form  an  Abelian  group,  we  may  choose 
a new  basis  for  the  group  from  among  products  of  the  oid.  Specifically, 


i EuRe  (j^f1  - fRe 


(9) 


I D PL_P 

where  f - j (p  ls  the  fricti  on  factor  commonly  used  in  captiiary 

Pv 

flow.  It  is  important  to  note  that  the  group,  fRe,  arises  only  when 
inertia  is  negligible.  Langlols  (12)  describes  an  alternative  definition 
of  P*.  appropriate  when  Inertia  is  negligible,  which  leads  directly  to 
the  group  we  call  fRe.  Finally,  the  end  correction  procedure  allows  us 
to  treat  fully  developed  flow  only,  so  we  may  neglect  L/D.  The  suspension 
flow  at  fixed  fiber  concentration  is  thus  described  by  a relationship  of 
the  form: 


F(fRe,  De)  - 0 

(10) 

fRe  ■ fRe  (De) 

(ID 

VII.  EFFECT  OF  OIL  VISCOSITY  ON  SUSPENSION  VISCOSITY 

Figure  8 shows  that  unique  curves  of  fRe  versus  De  are  indeed 
obtained,  one  curve  for  each  of  the  three  fiber  concentrations.  The 
symbols  in  Figure  8 represent  +6%  error  flags  (one  standard  deviation). 
We  believe  the  increased  scatter  in  the  30  vol .%  data  is  due  to  the  small 
change  in  the  viscosity  of  the  oil  which  occurred  during  mixing— -the 
0 vol.-  suspensions  were  not  mixed  and  the  15  voi.%  suspensions  were 
diluted  with  pure,  undegraded  oil,  so  the  change  in  viscosity  shows  up 
mostly  in  the  30  vol. 3:  data. 
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We  note  that 
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so  that  Figure  8 may  be  interpreted  as  a plot  of  reduced  apparent 
viscosity  versus  reduced  apparent  shear  rate.  The  curves  have  the 
general  form  normally  assumed  by  such  plots  for  polymeric  materials. 

The  only  theory  in  the  literature  which  predicts  quantitatively  the 
shear  rate  dependence  of  the  viscosity  of  concentrated  fiber  suspensions 
is  that  of  Ziegel  (4).  His  theory  does  not  fit  our  data. 


VIII.  EFFECT  OF  FIBER  CONCENTRATION  ON  SUSPENSION  VISCOSITY 

The  concentration  dependence  of  suspension  viscosity  is  normally 
presented  by  plotting  reduced  viscosity  against  concentration.  As  shown 
above,  reduced  apparent  viscosity  is  an  appropriate  variable  for  describing 
our  suspensions  but  it  must  be  expected  to  depend  on  Deborah  number.  The 
reduced  apparent  viscosi ty-versus-concentration  data  for  our  suspensions 
are  plotted  H Figure  9 for  various  Deborah  numbers.  The  reduced  apparent 
viscosity  is  less  than  one  at  h.gh  De  Because  it  is  referred  to  the  zero- 
shear-rate  viscosity  of  the  oil.  The  curves  resemble  those  reported  by 
Newman  and  Trementozzi  (13)  for  different  fiber  suspensions. 

As  can  be  seen  in  Figure  9»  the  dependence  of  suspension  viscosity 
on  concentration  is  highly  shear  rate  dependent.  > A qualitative  understanding 
of  the  strong  shear  rate  dependence  may  be  gained  from  the  dimensional 


n 


• • 


analysis  discussed  earlier.  It  is  obvious  from  Figure  8 that  the 
dimensionless  group  fRe  depends  on  De.  In  other  words,  the  silicone  oil 
is  not  behaving  as  a Newtonian  fluid  despite  that  fact  that  the  pure 
oil  is  Newtonian  over  most  of  the  flow  rate  range  considered.  This 
result  suggests  that  the  apparent  shear  rate  at  the  wall  Is  a poor  (low) 
measure  of  the  actual  shear  rates  occurring  In  the  suspension,  because 
the  increased  viscosity  of  a suspension  over  that  of  the  pure  fluid  Is 
associated  with  the  disturbed  velocity  profile  and  higher  shear  rates 
in  the  fluid  caused  by  the  presence  of  the  particle.  Presumably  these 
higher  shear  rates  lower  the  non-Newtonian  viscosity  of  the  fluid 
around  the  particle.  Hence,  the  resulting  suspension  viscosity  Is  lower 
than  would  be  expected  if  the  fluid  were  Newtonian.  This  effect  should  be, 
and  is,  more  pronounced  at  higher  fiber  concentrations  and  higher  flow 
rates. 

It  appears  in  principle  that  the  phenomenon  could  be  strong  enough 
that  the  apparent  viscosity  of  a suspension  might  lie  below  the  apparent 
viscosity  of  the  pure  oil  at  the  same  flow  rate.  Figure  8 suggests  this 
might  occur  in  our  suspensions  at  shear  rates  slightly  above  the  highest 
we  were  able  to  measure.  The  phenomenon  could  help  explain  the  data  of 
Thomas  and  Hagen  (14),  who  observed  suspension  viscosities  which  were 
lower  than  that  of  pure  resin  and  who  explained  the  results  as  due  to 
resin  degradation.  In  a recent  paper  Charrier  and  Rieger  (15)  found 
very  weak  dependence  of  viscosity  on  fiber  concentration  in  comnercial 
alass-f i I led  polypropylene,  and  they  have  invoked  a similar  mechanism  to 
explain  their  observations. 
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The  suspension  theories  of  Brodnyan  (3),  Ziegel  {k) , Hashin  (5) 
and  Nielsen  (6)  provide  equations  which  predict  the  dependence  of  reduced 
suspension  viscosity  on  fiber  concentration.  Zlegel's  theory  predicts 
linear  dependence  of  n on  $,  which,  as  Figure  8 shows,  is  not  exhibited 
by  our  suspensions.  Comparison  of  the  remaining  equations  with  our  data 
is  difficult  because  these  equations  assume  mono-d is perse  fibers  and 
Newtonian  oil  behavior.  We  chose  the  median  aspect  ratio  of  20:1  to 
represent  the  length  of  the  fibers.  The  maximum  allowed  volume  fraction 
for  these  fibers,  as  determined  from  centrifuging  experiments,  was  38  vol.fc. 
fhe  suspensions  appear  to  exhibit  a generalized  Einstein  coefficient,  k, 
of  11.5-  The  equations  using  these  values  are  plotted  in  Figure  10  and 
are  compared  with  the  data  taken  at  the  lowest  shear  rate,  because  the 
silicone  oil  is  behaving  most  like  a Newtonian  fluid  at  low  shear  rates. 
Brodnyan's  prediction  at  30  vol.9’  is  ten  orders  of  magnitude  too  high  and 
clearly  does  not  fit  the  data.  A similar  conclusion  was  reached  by 
Newman  and  Trementozzi  (13).  Hashin's  prediction  is  too  low,  but  the 
use  of  his  equation  with  our  data  is  questionable  because  our  suspensions 
dc  not  meet  his  assumptions  of  parallel,  infinitely  long  fibers.  Nielsen's 
equation  fits  the  data  within  experimental  uncertainty. 


IX.  DISCUSSION  OF  THE  DIMENSIONAL  ANALYSIS 

The  dimensional  analysis  presented  in  this  paper  is  not  rigorous. 

The  fact  that  it  makes  some  interesting  predictions  about  suspension  flow 
which  fit  our  data  suggests  that  a more  rigorous  dimensional  analysis 
might  be  worthwhile. 
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The  very  general  analysis  may  be  applied  to  viscometers  of  any 
geometry  and  to  particles  of  any  shape  In  suspensions  of  any  concentration. 
However,  the  group  fRe  a n /n  arises  only  when  Inertia  is  negligible, 
which  suggests  that  the  reduced  viscosity  Is  an  appropriate  parameter  for 
characterizing  suspensions  only  In  this  case.  The  analysis  also  predicts 
that  the  shear  rate  dependence  of  n /n  arises  from  De;  that  is,  from  the 
non-Newtonian  behavior  of  the  suspending  oil  and  not  from  changes  In 
particle  orientation  as  is  commonly  asserted. 

However,  two  problems  exist  with  the  analysis  as  developed. 

First,  the  analysis  Is  based  on  an  assumption  about  the  initial 

condition  of  the  suspension  in  the  barrel  but  the  results  of  the 

ana • ys is  are  applied  to  fully  developed  flow.  This  slight  of  hand  only 

presents  significant  problems  when  we  wish  to  neglect  inertia.  A good 

argument  c«m  be  made  for  neglecting  inertia  In  fully  developed  flow; 

what  is  needed  is  a clear  statement  of  what  is  meant  by  neglecting  inertia 

during  start-up  in  the  entry  region  of  the  viscometer. 

A more  rigorous  statement  of  the  basic  assumption  of  statistically 
equivalent  distributions  of  initial  fiber  positions  Is  also  needed.  Some 
insight  into  the  problem  may  be  gained  by  considering  the  ensemble  of  all 
possible  distributions  a'  represented  in  Figure  II  for  two  suspensions, 

A and  B,  which  differ  only  in  oil  viscosity.  Both  ensembles  are  identical 
because  the  initial  fiber  positions  which  are  possible  do  not  depend  on 
the  oil  viscosity.  Our  data  tells  us  that  each  ensemble  contains  a large 
subset  comprising  the  initial  distributions,  the  members  of  which  are 
experimentally  indistinguishable.  These  initial  distributions  are 


■ '...I. 


Hi 

contained  in  the  curves  os  Figure  1). 


Typical  extrusion  force  measurements 


on  suspensions  A and  B will  be  made  on  some  distributions  which  are  members 

j 

of  the  subset--call  them  Aj  and  Bj . The  existence  of  these  subsets  are 
fundamental  to  reproducible  viscosity  measurements  on  suspensions  v/lthout 
explicit  control  over  particle  orientation.  The  exact  nature  of  the 
subsets  presumably  depends  on  the  accuracy  of  the  experiments.  We  assume 
that  the  two  subsets  have  at  least  one  common  member--distrlbutlon  0 in 
Figure  II.  This  assumption  is  reasonable  in  view  of  the  large  number  of 
distributions  in  each  subset.  The  results  of  measurements  on  Aj  and  Bj 


are  experimentally  indistinguishable  from  A and  B , and  the  latter  mav 

o o ' 

be  rigorously  compared  by  dimensional  analysis.  Thus  we  have  some 
confidence  that  wc  can  treat  members  such  as  Aj  and  Bj  by  the  same 
analysis. 


1 

1 
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X.  CONCLUSIONS 

This  paper  described  an  experimental  technique  by  which  several 
suspensions  with  the  same  fiber  length  distribution  but  different  oil 
viscosities  and  fiber  concentration  were  mixed.  The  following  conclusions 
apply  to  those  suspensions  in  the  range  in  which  the  measurements  were  made. 


1 . 


Wall  effects  and  other  special  suspension  effects  were  not 
detected . 


2.  Suspension  viscosity  was  not  dependent  on  capillary  entrance 
angle  or  capillary  diameter. 
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3.  Suspension  viscosity  increased  with  oil  viscosity  and  generally 
decreased  with  flow  rate.  The  data  fit  a d.' mens  ion  less  plot  of 
the  form 

F(fRe,  De)  - 0 

which  is  essentially  a plot  of  reduced  apparent  viscosity  versus 
reduced  apparent  shear  rate. 

4.  Suspension  viscosity  increased  with  fiber  concentration,  and  the 
magnitude  of  the  increase  was  highly  shear  rate  dependent.  One 
explanation  of  the  shear  rate  dependence  is  the  non-Newtonian 
behavior  of  the  silicone  oil. 
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Table  ) 

Fiber  Length  Distributions  of  the  Suspensions 
(Uncertainties  are  two  standard  deviations) 


Oil 

Viscosity 

(P) 


Concentration  Mean  Fiber 

(y°l '%)  Length 

Oog^  Inch) 


Standard 
Deviation 
(log^  Inch) 


100 

100 

600 

600 

1000 

1000 


15 

30 

15 

30 

15 

30 


"5-38  ♦ 0.08 
"5.1*9  + 0.08 
"5.30  + 0.08 
-5.26  + 0.08 
-5.21*  + o.08 
-5.27  +•  0.08 


0.71  ♦ 0.07 
0.71  + 0.07 
0.67  ♦ 0.07 
0.66  ♦ 0.07 
0.68  ♦ o.(J7 
0.66  ♦ C.07 


I'i 

Table  2 

Suspension  Parameters  which  Determine  Suspension  Viscosity 

Material  Properties: 

Oil  Rheology: 

Viscosity  - 100,  600,  1000  P 
Time  dependent  and  elastic  effects  - negligible 
Fiber  Properties: 

Length  - same  distribution  in  all  suspensions 
Diameter  - 0.0005" 

Youngs  modulus  - 10^  psi 
Poisson  ratio  - 0,2 

Suspension  Properties: 

Concentration  - 0,  15,  30  vol.fc 
Disperson  and  Orientation  - random 
Third  phases  - none 
Surface  effects  - negligible 


II 


Entrance 

Angle 

(°) 


180 

90 


Nominal 

Viscosity 

(P) 

100 

600 


Table  3 

Effect  of  Viscometer  Entrance  Angle 
(Uncertainties  are  One  Standard  Deviation) 


Shear 

Rate 

(sec"') 


527 

527 


Entrance 
Pressure 
Losses  2 
(dynes/cm  ) 


6.3  + 1.3  x 106 

7.3  + 0.6  x 106 


Viscosity 

(P) 


960  ♦ 36 
940  + 20 


Table  4 

Carreau  Model  Constants  for  the  Silicone  Oils 


o 

o 

(P) 

(sec"’ ) 

100 

2300 

0.70 

600 

400 

0.67 

1000 

240 

0.67 

1000 


CAPTIONS  FOR  FIGURES 


OP 


1.  Histogram  rf  the  fiber  lengths  in  a typical  suspension. 

2.  Shear  stress-shear  rate  behavior  of  three  suspensions  of  fibers  in 
100  P oil.  The  symbols  (Q-l/4"  cap.;/^-l/8"  cap.)  represent  +6$ 
error  flags. 

3.  Shear  stress-shear  rate  behavior  of  three  suspensions  of  fibers  in 
600  P oil.  The  symbols  (0"'A"  cap.;A-l/8"  cap.)  represent  +6$ 
error  flags. 

4.  Shear  stress-shear  rate  behavior  of  three  suspensions  of  fibers  in 
1000  P oil.  The  symbols  (O  "1A"  cap.;  A-1/8"  cap.)  represent 
+6?;  error  flags. 

5.  Viscosity-shear  rate  behavior  of  three  suspensions  of  fibers  in 
100  P oil.  The  symbols  (0“lA"  cap.;  ^-1/8"  cap.)  represent 
+6%  error  flags. 

6.  Viscosity-shear  rate  behavior  of  three  suspensions  of  fibers  In 
600  P oil.  The  symbols  (0"'A"  cap.;  ^-1/8"  cap.)  represent 
+6%  error  flags. 

7.  Viscosity-shear  rate  behavior  of  three  suspensions  of  fibers  in 
1000  P oil.  The  symbols  (Q-l/4"  cap.;  ^-1/8"  cap.)  represent 
+6%  error  flags. 

8.  Dimensionless  plot  of  fRe  versus  De  at  constant  n.  The  symbol 
size  represents  +69;  error  flags. 


9.  Reduced  viscosity-concentration  dependence  of  all  suspensions  at 
several  Deborah  numbers.  The  symbols  have  no  significance  as 
error  flags. 


10.  Comparison  of  viscosity-concentration  dependence  at  lowest  Deborah 
number  with  literature  predictions. 


T 11.  Representations  of  the  ensembles  of  all  possible  initial  fiber 

I position  distributions  for  two  suspensions  with  different  oil 

viscosities.  Each  small  square  represents  a particular  distribution. 

JThe  distributions  inside  the  closed  curve  are  experimentally 

indist i ngu i shable . 
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?lot  of  the  dimensionless  groups  of  the  system  et  constant 
represents  : A error  flags. 


concentration  dependence  at  lowest  Deborah  number 


Representations  of  the  ensembles  of  all  possible 
Initial  fiber  position  distributions  for  two  suspensions 
with  different  oil  viscosities.  Each  small  square 
represents  a particular  distribution.  The  distributions 
inside  the  closed  curve  are  experimentally  indistinguishable. 


